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ABSTRACT 

A  study  of  single-span  initially  bent  columns  is  made  and  extended 
to  include  two-span  columns.   The  axial  force  in  these  columns  results 
from  heating  (or  an  equivalent)  and  from  bowing.   A  snap-through 
phenomenon  is  considered  in  the  analysis  of  the  single-span  and 
observations  are  made  concerning  snap-through  in  the  two-span  column. 
A  digital  computer  program  is  developed  to  carry  out  the  analysis  and 
a  description  of  the  program  enables  its  use  by  anyone  familiar  with 
the  FORTRAN  language. 
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I.   INTRODUCTION 

A.   GENERAL  REMARKS 

The  subject  of  initially  bent  columns  and  shallow  arches  is 
treated  extensively  in  the  literature.   Hsu  [Ref.  6]  studied  the 
dynamic  stability  of  shallow  arches  against  snap- through ;  Hoff  and 
Bruce  [Ref.  5]  analyzed  dynamically  the  snap-through  of  shallow 
arches;  Fung  and  Kaplan  [Ref.  2]  studied  the  buckling  of  low  arches 
or  curved  beams  of  small  curvature;  Timoshenko  [Ref.  8]  has  also 
discussed  the  problem;  Huddleston  [Ref.  7]  studied  the  effect  of 
axial  strain  in  elastic  columns;  Gjelsvik  and  Bodner  [Ref.  3] 
approached  the  problem  on  the  basis  of  energy  criteria.   There 
have  been  other  similar  studies  made.   However,  these  treatments 
of  the  problem  have  all  involved  only  a  single  span.   Brock  [Ref.  1] 
studied  the  single  .span  using  Fourier  analysis  in  such  a  way  as  to 
provide  for  assembly  into  a  multispan  continuous  column,  and  gave 
conditions  for  the  solution  of  multispan  cases.   However,  he  did 
not  actually  arrive  at  numerical  solution  in  any  particular  case, 
only  suggesting  how  the  computations  might  be  carried  out.   There- 
fore, it  v/as  the  objective  of  this  thesis  not  only  to  study  the 
single-span  column  in  greater  depth,  but  to  develop  its  analysis 
in  a  manner  such  that  it  could  most  easily  be  extended  to  a  multi- 
span  problem,  and  actually  to  make  this  extension  to  the  case  of  a 
two-span  column. 


B.  SCOPE  OF  WORK 

The  theory  of  the  single-span  problem  is  first  developed  for  a 
column  which  is  initially  nearly  straight.   This  development  takes 
into  consideration  an  axial  strain  induced  in  the  column  and  also 
allows  for  the  external  application  of  end  moments.   It  is  assumed 
that  the  configuration,  with  and  without  load,  lies  in  a  single 
plane.   The  theory  is  then  extended  to  include  the  two-span  column. 
A  digital  computer  program  is  developed  to  carry  out  the  analysis, 
and  a  description  of  the  program  enables  its  use  by  anyone  who  is 
familiar  with  the  FORTRAN  language. 

Analysis  beyond  the  two-span  problem  is  not  carried  out  in  this 
thesis,  although  it  is  envisioned  that  multi-span  problems  can  be 
investigated  in  much  the  same  manner  as  described  herein. 

C.  NOTATION 

A  -  cross  sectional  area 

a  -  Fourier  coefficients  for  initial  configuration 

b  -  Fourier  coefficients  for  final  configuration 

E  -  Young's  modulus  of  elasticity 

I  -  moment  of  inertia  of  column 

L  -  length  of  span 

M  -  left  end  moment 

N  -  right  end  moment 

n  -  term  index  in  Fourier  analysis 

P  -  axial  compressive  load 

P  -  Euler's  critical  buckling  load 
c  ^ 

Q  -  left  end  reaction  force 
R  -  right  end  reaction  force 


W  -  lateral  load  in  two-span  column 

X  -  distance  from  left  end  of  span 

V  -  initial  deflection 
•^o 

y  -  final  deflection 

e  -  unit  strain 

e  -  thermal  unit  strain 

K  -   curvature 

X   -  excess  of  arc  length  over  chord  length  in  initial 
configuration 

A   -  excess  of  arc  length  over  chord  length  in  final 
configuration 

71   -  3.141592  ... 

6   -  end  slope  of  span 


II.   DEVELOPMENT  OF  THEORY 

A.   SINGLE  SPAN 

The  study  of  the  single  span  involved  an  initially  bent  column 
loaded  as  shown  below: 


Figure  1.   Loading  on  single  span  column 


The  end  supports  are  frictionless  hinges  which  do  not  move  from 
their  original  positions.   The  moment  loading  shovm  is  applied  by 
external  agencies.   The  initial,  no-load  configuration  of  the  column 
is  presumed  to  be  given  by  the  equation 


V  =  E  a  sin 
■^o      n 


nrrx 

L 


(1) 


where  here  and  later  the  summation  is  viith   respect  to  n  and  goes  from 
n=l  to  n=  •» .   For  computational  purposes  the  upper  limit  is  taken  to 
be  a  finite  number,  which  will  be  discussed  in  a  later  chapter. 

Next  suppose  that  the  column  is  uniformly  heated  so  as  to  provide 
thermal  strain.   (It  should  be  noted  that  an  equivalent  situation 
would  result  from  moving  the  ends  of  the  column  tov/ard  each  other. 
However,  throughout,  we  will  speak  in  terms  of  heating  the  column 
rather  than  displacing  its  ends.)   The  new  configuration  can  be 
presumed  to  be  represented  by  the  equation 

(2) 


y   =^\   sin(l^). 
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The  immediate  purpose  is  to  find  the  b  's  and  thereby  the  new 
configuration. 

There  is  an  axial  compressive  force  P  set  up  in  the  column  such 
that  the  unit  strain  is 

e  =  E^  -  P/AE  (3) 

where  ^   is  the  thermal  strain,  A  is  the  cross  section  area,  and  E 

is  Young's  modulus  of  elasticity.   (Throughout,  we  presume  that  Hooke'; 

law  is  satisfied. ) 

The  excess  of  arc  length  over  chord  length  can  be  adequately 
represented  for  the  initial  deflection  by  the  approximate  relation 

o    2   /   dx       4L        n 

0 
Likewise  for  the  new  configuration 


X  =  ^  5:  n^b  2.  (5) 

4L      n 


Thus,  the  change  in  length  due  to  heating  is 

^  -^   =  (ttV   n^(b  2  -  a  2)  =  Le^  -  PL/AE.  (6) 

o   V^L/      n     n       T 

The  bending  behavior  of  the  column  must  also  be  accounted  for.   By 

static  equilibrium 

Q  =  -R  r  (N-M)/L.  (7) 

The  bending  moment  at  any  point  is  given  by 

B.M.  =  QX  +  M  -  Py,  (8) 

and  this  accounts  for  the  change  in  curvature 

A^  d2y      ,  .  2 

d  y     -^o   '-^ 


Ak  -  —Y  ~  -7—2  - 
dx    dx 


(f)  ^"^Cb„-a„)sin(nSi).  ,3, 


According  to  the  elementary  formula, 

B.M.  =  EI  Ak.  (10) 

Thus,  with  a  little  manipulation,  we  find 

^    [n^P  (b   -  a  )  -  Pb  ]  sin  (^  }+  Qx  +  M  =  0 ,  (11) 

cnn      n      ^Ly  ' 

where  we  have  introduced  the  notation 

P   =   ■^^EI/L^  .  (12) 

c 

This  is  the  Euler  critical  load  which  appears  naturally  in  this 

calculation.   To  deal  successfully  with  Eq.  (11)  it  is  useful  to 

express  Qx  and  M  as  Fourier  series  in  the  range  0  _<  x  <_  L. 

Qx  3  (2QL/^)  ^(-l)""^^  (1/n)  sin  h~]  (13) 

M  =    (2M/TT)  i:[l  -  (-i)""]  (1/n)  sin  (^  ]   .  (14) 

Thus, 


0  =  l<i-l)^^ 


1  f  2QL  \  ^  ^i_(.i)n  /  2M  \  ^  ^2p  (^  _,  )_p,  1  3^,  M,      (15) 
\    m^  j  \"/      cnn    n/     \    ^    / 


which  will  be  satisfied  if  each  of  the  terms  in  braces  vanishes. 
Accordingly,  we  set 


^^  "\V^->"(H^^<-'"-^(l^W 


P) 
c 


n=  1,2,...   (16) 


However  this  equation  still  involves  the  unknown,  P.   Equation  (5)  can 
be  written  in  the  form 

P  =  AEe^  -Ti^AE/UL^  Z    n^  (b  ^  -  a  2).  (17) 

T  n    n 

Rewriting  Eq.  (17)  once  more  in  the  form 

f(P)  =  0,  (18) 

where 

10 


f(P)  =  PL/AE  +  (TT^/i+DSn^  (b  2  _  a  2)  -  Le^,  (19) 

n     n       T 

we  seek  the  zero  (or  zeros)  of  the  nonlinear  function  f(P).   Because 
of  nonlinearity ,  computational  difficulties  were  encountered  in  find- 
ing the  significant  root.   Of  several  methods  attempted,  the  most 
suitable  was  the  Newton- Raphson  iteration  scheme.   A  study  of  the 
curve  f(P)  revealed  that  the  method  would  usually  converge  on  the 
first  root  if  the  iteration  was  started  at  a  value  just  slightly 

smaller  than  P  . 
c 

The  value  of  P  thus  being  known,  Eqs.  (16)  can  then  be  solved  for 

the  coefficients  b  . 

n 

Figure  2  shows  the  initial  and  final  deflection  curves  of  a  column 
computed  using  this  analysis.   The  final  deflection  shows  the  column 
after  having  undergone  thermal  strain.   The  end  moments  were  zero  for 
this  simple  case. 

Further  analysis  of  the  single  span  problem  will  be  given  in 
Section  III  of  the  thesis. 

B.   TWO  SPAN 

We  now  consider  two  such  spans,  the  column  being  joined  at  the 
intermediate  support  so  as  to  maintain  continuity  of  slope  and  moment 
as  well  as  of  deflection.   In  general,  different  axial  loads  will  be 
acting  in  the  two  spans.   As  before,  the  problem  is  to  determine  the 
final  configuration  after  heating,  the  initial  configuration  having 
been  specified. 

The  obvious  method  of  analysis  is,  in  effect,  to  consider  two 
adjacent  single  spans,  each  of  which  is  as  shown  in  Fig.  1  and  is 
analyzed  as  in  the  preceeding  subsection.   The  adjoined  spans  are  as 
shown  in  Fig.  3.   Separating  them  and  considering  their  mutual 
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interaction  and  the  continuity  at  their  junction,  we  have  (see 
Fig.  4) 

^1  =  S' 


(20) 


Rl    L2  • 


(21) 


If  the   coiranon   value  N      =   M     can  be   determined,    the  problem  is 

solved  in  the  sense  that  each  span  can  be  now  analyzed  separately.   As 

a  matter  of  fact,  in  the  solution  for  N   =  M  ,  we  will  actually  have 

determined  P,  and  P.  and  the  Fourier  coefficients  b   for  each  span. 
12  n 


R.*0. 


Figure  3.   Loading  on  two-span  column 


Figure  4.   Separation  into  two  single-span  columns 


The  left  and  right  end  slopes  of  each  span  may  be  found  by 

evaluating  dy/dx  at  x  =  0  and  x  =  L,  respectively. 

nfTb 

nTTb 

'r  =   ^(-1)"  -T^  • 


(22) 
(23) 
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A  method  of  solution  for  the  two  span  problem  is  outlined  as 
follows.   First  assume  a  common  value  for  N  and  M  .   (This  will  be 
called  the  junction  moment.)   Then  analyze  each  span  separately,  using 
the  analysis  of  the  preceding  section,  and  determine  the  slopes  at 
the  junction.   Compare  the  two  slopes  and  find  the  difference  between 
them,  which  will  be  called  error.   Next  repeat  the  above  for  another 
value  of  moment  and  find  the  error.   A  plot  of  error  versus  moment 
will  indicate  the  value  of  junction  moment  that  will  cause  the  error 
(in  slope)  to  be  zero.   At  this  point,  correct  values  for  the  Fourier 
coefficients  b  will  be  known  for  each  span  so  that  the  final  config- 
uration can  then  be  computed  and  plotted  if  required. 

If  the  plot  of  slope  error  versus  junction  moment  is  nearly 
linear,  a  slope  intercept  method  of  iteration  may  be  used  to  determine 
where  the  error  goes  to  zero.   This  method  simply  takes  two  points, 
gets  the  equation  of  the  line  passing  through  these  points,  and  then 
solves  for  the  zero  of  this  linear  approximation.   This  is  repeated, 
using  a  new  point  and  one  of  the  old  points,  until  the  desired 
convergence  is  reached.   If  the  plot  of  slope  error  versus  junction 
moment  is  significantly  nonlinear  (this  occurs  if  the  axial  compressive 
load  is  a  large  fraction  of  the  Euler  load),  then  the  slope  intercept 
method  is  no  longer  satisfactory  and  a  method  such  as  successive 
decimation  must  be  used.   In  this  context,  decimation  means  recalcu- 
lating over  an  interval  of  interest,  using  increments  one  tenth  as 
large  as  previously.   The  program  decimates  four  times  so  that  the 
final  increment  is  .0001  times  the  initial  increment.   This  method  is 
well  known  and  will  not  be  described  here.   The  accuracy  of  results 
obtained  by  using  this  method  depends  upon  the  size  of  the  increment 
used  in  stepping.   The  smaller  the  increment,  the  greater  the  accuracy. 
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The  question  here  is  where  to  start  the  iteration.   If  a  good 
approximate  figure  is  known,  then  the  increments  may  be  small. 
If  not,  the  increments  must  be  larger.   Unfortunately  there  is  no 
quick  way  of  determining  a  good  starting  approximation. 

C.   INITIAL  DEFLECTION 

For  the  single  span  problem,  any  initial  configuration  may  be 
specified  subject  only  to  the  condition  y(0)=y(L)=0.   If  the 
configuration  is  specified  in  the  form  of  a  polynomial  in  x,  the 
axial  distance  from  the  left  end,  the  procedure  to  be  described 
later  quickly  determines  as  many  of  the  Fourier  coefficients  a   as 
may  be  desired. 

However,  for  the  two  span  problem  the  necessity  for  continuity 
places  extra  requirements  on  the  initial  deflection.   Not  only  must 
the  deflection  be  zero  at  each  end  of  the  spans,  but  the  slope  and 
curvature  must  be  continuous  at  the  junction.   Accordingly,  such  a 
two  span  configuration  is  obtained  as  follows.   An  arbitrary  config- 
uration (conveniently  in  the  form  of  a  polynomial)  is  assumed  which 
has  zero  displacement  at  the  left  of  span  one  and  at  the  right  of  span 
two.   (see  Fig.  5).   Then,  regarding  this  as  a  simply  supported  beam, 
sufficient  lateral  load  is  applied  at  the  intermediate  support  to 
cause  the  displacement  to  vanish  at  this  point.   The  corresponding 
deflection  curve  is  calculated  and  its  superposition  with  the 
originally  assumed  configuration  gives  a  result  which  now  has  zero 
deflection  at  each  of  the  three  support  points.   In  general,  the 
column  will  thus  start  out  with  non-zero  bending  moment  except  at 
the  extreme  ends. 


15 


y  =  f(x) 


W  (to  be  applied) 


Figure  5.   Original  specified  configuration 


Figure  5  illustrates  the  originally  specified  configuration 
before  the  intermediate  constraint  is  enforced.  It  is  not  dif- 
ficult to  show  that  the  necessary  lateral  load  is  given  by 


Force  =  W  = 


SCL^+L^)  EIA 
2 2 

4S 


downward. 


(24) 


It  can  also  easily  be  shown  that  the  deflection  of  a  beam  due  to  a 
downvjard  force  W  at  point  B  is  represented  by 


6EI(L  +L  ) 

4_ 


,         (L.+L  ) 


(L^+L^)^  x-x  -  L^^x  + 


<x  -  L  >3 


where  the  <  >  symbol  has  the  meaning 


(25) 


<  X  -  L^>' 


if  x  <  L. 


(x-L^)^   if  X  1  L^ 


Thus,  by  superposition,  the  new  configuration  is 


WL, 


y  =  ^^^^  -  &etu-^ 


(L^+L^) 
(L^+L^)^  X  -  x^  -  L^^x  + 


<x  - 


h' 


(25) 


For  the  first  span  the  expression  for  the  configuration  is 


where  0  <  x  <  L  . 


(L^+L^)^  -  x2  -  L^^ 


C27) 


15 


For  the  second  span 

WL^Cx-'HL^) 
y  =   f(x:':+Lj_)  -  s^j(L.^^L^) 


(L  +L  ) 


(28) 


where  X"  =  x-L^  is  the  distance  along  span  two  from  its  ovm  left  end. 
Now  the  Fourier  coefficients  can  be  found  for  each  span.  If  f(x) 
is  a  polynomial,  the  program  described  in  Appendix  A  makes  the  deter- 
mination. Thus  we  have  a  deflection  satisfying  the  necessary  initial 
conditions  for  the  two  span  problem. 
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III.   DISCUSSION 

A.   ANALYSIS  OF  THE  PROBLEM 

In  analysis  of  the  single-span  problem,  the  effects  of  the 
application  of  end  moments  were  investigated.   This  means  that 
end  moments  vjere  externally  applied  at  one  or  both  ends  of  the 
column  with  or  without  the  presence  of  axial  strain.   It  was 
found  that  a  snap-through  phenomenon  could  be  observed  using 
this  analysis.   That  is,  by  increasing  one  or  both  end  moments, 
the  configuration  could  be  seen  to  "jump"  from  one  configuration 
to  a  distinctly  different  configuration  with  only  infinitesimal 
variation  in  moment.   Figure  6  illustrates  this  phenomenon  by 
showing  a  column  that  is  first  heated  up,  and  then  subject  to  a 
left  end  moment.   Curve  1  represents  the  initial  configuration 
of  the  column.   Curve  2  is  the  configuration  of  the  column  after 
undergoing  thermal  strain.   The  remaining  curves  show  the  config- 
uration under  the  action  of  several  different  values  of  end  moment. 
It  is  seen  that  the  first  few  1000  inch  pound  increments  of  applied 
moment  cause  only  slight  changes  in  configuration,  but  that  the 
change  from  6000  to  7000  causes  a  relatively  large  change  which 
has  the  nature  of  a  jump,  or  snap-through  phenomenon,  following 
which  addition  of  increments  of  applied  moment  again  result  in 
only  minor  adjustments  of  the  configuration. 

Figure  7  further  illustrates  the  phenomenon  by  showing  a  plot 
of  centerline  deflection  versus  time,  where  time  is  used  in  the 
sense  that  first  thermal  strain  is  increased  uniformly,  with  no 
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end  inoment ,   and  then  thermal  strain  is  held  constant  and  the  end 
moment   is  increased  at  a  uniform  rate. 

It  may  be  noted  that  a  similar  analysis  could  be  made  for  a 
column  under  lateral  loading  instead  of  end  moment  loading. 

It  will  be  of  interest  to  explain  why  the  phenomenon  of  snap- 
through  is  considered  to  be  of  great  importance.   In  a  piping  system 
which  is  heated  up,  considerable  energy  is  stored  in  the  form  of 
strain  energy.   If  snap-through  or  a  similar  phenomenon  could  take 
place,  some  of  this  energy  might  be  suddenly  released  and  have  a 
catastrophic  effect  on  the  piping  system  or  adjacent  installations. 
In  the  case  of  a  single-span  straight-line  system,  if  the  ends 
should  be  pinned,  so  that  no  moments  may  be  applied,  it  is  clear 
that  snap-through  cannot  take  place.   If  the  ends  are  fixed  to 
heavy  equipment,  such  as  pumps  or  vessels,  the  restraining  influence 
of  these  items  would  indeed  result  in  end  moments  being  applied  to 
the  pipe,  but  it  is  difficult  to  envision  these  moments  actually 
resulting  in  snap-through. 

However,  in  the  case  of  multi-span  straight-line  configurations, 
at  the  intermediate  supports  adjacent  spans  exert  moments  on  each 
other,  and  it  is  not  immediately  clear  that  as  the  system  is  heated, 
perhaps  nonuniformly ,  these  intermediate  moments  might  not  increase 
in  such  a  way  as  to  cause  snap-through.   The  intuitive  feeling  that 
this  would  not  take  place  was  insufficient  basis  for  asserting  that 
it  could  not  happen.   Accordingly,  considerable  attention  was  focused 
on  this  possibility.   It  was  only  after  studying  many  different  initial 
configurations  and  thermal  loading  sequences  that  a  reasonable  under- 
standing was  achieved  of  the  restraining  influence  which  one  span  has 
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upon  an  adjacent  span  and  of  how  this  serves  to  prevent  the  phenomenon 
of  snap-through,  even  if  one  end  of  a  multi-span  configuration  is 
subjected  to  increasing  external  moment. 

Clearly,  the  simplest  multi-span  configuration  is  that  having 
two-spans.   Accordingly,  this  case  was  studied  in  detail  and  computer 
programs  were  written  for  its  analysis.   One  of  the  cases  which  was 
studied  is  illustrated  in  Fig.  8.   In  this  case  a  negative  moment 
of  increasing  magnitude  was  applied  at  the  extreme  left  causing  the 
column  to  go  through  the  succession  of  configurations  shown.   (Inci- 
dentally, although  initially  there  was  axial  compression  in  the  column 
due  to  initial  heating,  by  the  time  the  last  configuration  shown  was 
obtained,  the  axial  load  in  the  left  span  had  become  tensile.) 

This  example  provided  the  insight  to  see  why  snap-through,  at 
least  in  any  sense  similar  to  that  for  the  single-span  column,  is 
not  possible  for  two-span  columns.   The  original,  incorrect,  view 
was  that  moment  loading  at  the  left  end  of  a  two-span  column  was, 
in  effect,  simply  a  way  of  applying  moment  to  the  left  end  of  the 
right-hand  span;  increasing  the  leftmost  moment  sufficiently  V70uld 
eventually  apply  a  very  large  moment,  at  the  junction  point,  to  the 
left  end  of  the  right  span.   This  is  true  enough,  but  the  continuity 
between  left  and  right  spans  does  not  permit  the  latter  to  snap- 
through  as  would  be  the  case  if  it  were  the  only  span  and  the 
junction  moment  were  directly  applied  at  its  left  end.   In  effect, 
the  left  span  acts  as  a  stabilizing  influence  to  prevent  snap-through 
of  the  right  span. 

The  situation  may  be  visualized  more  clearly  by  imagining  that 
the  continuity  between  the  two  spans  is  "one-way"  in  the  sense  of 
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the  mechanism  shown  in  Fig.  9.   If  the  two  spans  were  joined  in  this 
manner,  snap- through  of  the  right  span  would  indeed  take  place  when 
the  junction  moment  reached  a  critical  value.   However,  with  the 
true  continuity  of  a  continuous  column  such  behavior  can  not  take 
place. 


M 


WITH    THIS    PIN    REMOVED,  AS 
SHOWN.  CONTINUITY    IS 
"ONE-WAY" 


RIGHT 
SPAN 


Figure  9.   "One-way"  continuity  of  two-span  junction 

Attempts  to  create  a  situation  in  v/hich  snap-through  takes  place 
in  a  two-span  column  have  proved  unsuccessful.   Study  of  the  behavior 
in  these  cases  has  added  strength  to  the  conviction  that  such  a 
phenomenon  cannot  take  place,  even  under  the  action  of  increasing 
externally  applied  moments  at  the  ends  or  at  the  intermediate  supports. 
However,  this  growing  conviction  cannot  be  understood  as  a  convincing 
demonstration  that  the  phenomenon  cannot  indeed  take  place,  perhaps 
only  for  certain  unfortunate  combinations  of  physical  parameters. 
Accordingly,  further  study  of  snap-through  in  multi-span  columns  is 
suggested  as  an  appropriate  topic  for  further  investigation. 
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B.   FOURIER  REPRESENTATION  OF  THE  CONFIGURATIONS 

A  question  arises  of  the  ability  of  Fourier  sine  series,  as 
employed  herein,  to  represent  non-zero  curvature  at  the  ends  of 
spans  where  each  sine  term  vanishes.   Certainly,  the  value  of 
curvature  given  by  such  a  series  must  be  zero  when  x=0  and  x=L. 
If  we  have  initial  deflections  which  have  non-zero  curvature  at 
the  ends  of  the  span  (or  spans),  our  Fourier  series  must  be  able 
to  represent  this  adequately.   Also,  if  v;e  have  non-zero  moment 
at  these  points  (due  to  the  application  of  external  moments  or  due 
to  enforcing  an  intermediate  constraint)  this  corresponds  to  non- 
zero curvature  change,  which  must  be  adequately  represented  by  our 
Fourier  series. 

Accordingly,  actual  plots  of  curvature,  as  determined  from 
the  Fourier  series ,  were  made  in  several  cases  where  the  value  of 
curvature  could  be  determined  by  other  means.   Figures  10,  11,  and 
12  show  a  comparison  of  curvature  in  a  case  where  the  initial 
configuration  (in  the  sense  of  Fig.  5)  had  uniform  negative 
curvature  and  where  the  intermediate  constraint  produced  a  piece- 
wise  linear  moment  change  (i.e.,  curvature)  that  v/as  zero  at  each 
end  and  a  maximum  at  the  constraint.   The  final  curvature  is  the 
superposition  of  these,  consisting  of  two  straight  segments  rising 
from  -.001332  to  .000665  and  back  down  to  -.001332.   (In  Figs.  10, 
11,  and  12  these  straight  line  segments  show  a  very  small  waviness 
due  to  the  method  of  plotting.)   These  figures  also  show  the  Fourier 
series  representation  of  curvature,  using  25,  50,  and  100  terms, 
respectively.   It  is  seen  that  although  indeed  the  Fourier  repre- 
sentation of  curvature  does  vanish  at  the  constraint  points,  this 
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is  a  very  local  phenomenon  and  elsewhere  the  Fourier  series  provides 
a  good  representation  of  the  correct  values. 

Another  example  is  shown  in  Fig.  13  where  values  of  curvature 
as  represented  by  Fourier  series  are  plotted  for  different  values 
of  left  end  moment  in  a  case  of  snap- through.   Fairing  straight 
lines  through  the  left  ends  of  the  curves  gives  intersections  with 
the  left  vertical  axis  which  correspond  to  the  actual  moments  applied 
at  the  left. 

Probably  the  most  convincing  verification  of  the  adequacy  of  the 
Fourier  representation  is  the  relative  insensitivity  of  significant 
results  (in  a  particular  problem)  to  the  number  of  Fourier  terms 
employed.   In  a  particular  case  of  snap-through,  Table  I  below  shows 
that  snap-through  occurred  in  the  same  moment  interval  for  all  values 
of  n  down  to  as  lovi   as  3.   It  is  interesting  to  note  in  Table  II  that, 
even  for  n=l,  the  axial  load  obtained  is  in  very  good  agreement  with 
the  load  obtained  using  100  terms. 


Table  I.   Values  of  end  moment  between  which  snap-through 
occurred,  for  different  values  of  n. 


Moment  (in  .lb.) 

n=100 

n=50 

n=25 

n=10 

n=5 

n=3 

n=2 

n=l 

60450 

• 

• 

• 

• 

• 

• 

60460 

• 

9 

60470 

Note:   Dots  indicate  the  moment  interval  in  which  snap-through 
occurred.   The  tabulation  is  for  a  particular  problem 
described  on  page  31. 
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Table  II.   Value  of  axial  load  In  pounds  for  different 
values  of  n  and  M. 


Moment  (in. lb.) 

60450 

60460 

60470 

n=100 

38489.000514 

38491.906680 

38489.899525 

n=50 

38489.000515 

38491.906680 

38489.899526 

n=25 

38489.000517 

38491.906680 

38489.899527 

n=:10 

38489.000550 

38491.906674 

38489.899548 

n=5 

38489.000764 

38491.906634 

38489.899684 

n=3 

38489.009402 

38491.905053 

38489.905163 

n=2 

38489.011286 

38491.009623 

38489.906358 

n=l 

38489.011298 

38491.009624 

38489.906349 

The  tabulation  is  for  a  particular  problem  described  below. 

Due  to  computational  difficulties  in  the  vicinity  of  snap- 
through,  values  of  moment  taken  at  smaller  than  10  inch-pound  in- 
crements provided  somewhat  erratic  results. 

The  results  tabulated  in  Tables  I  and  II  are  for  a  steel  pipe 
whose  cross  sectional  area  is  2.945  in^.   The  length  is  200  inches, 
the  moment  of  inertia  of  the  cross  section  is  5.2  in  ,  and  Young's 
modulus  of  elasticity  is  30  x  10°psi.   Moment  is  applied  at  only 
one  end  of  the  span.   The  initial  configuration  was  similar  to  that 
shown  in  Fig.  6. 
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IV.   ANALYSIS  or  A  TYPICAL  TWO-SPAN  PROBLEM 

Using  the  program  described  in  Appendix  A,  configurations 
resulting  froin  the  application  of  many  different  combinations  of 
end  moments  and  strain  on  a  two-span  column  can  be  obtained.   The 
following  cases  illustrate  several  of  the  ways  in  which  the  program 

may  be  utilized.   These  examples  involve  a  two-span  steel  pipe  with 

.2  .       . 

cross  sectional  area  of  2.945  m  .   The  left  span  is  100  inches 

long  and  the  right  span  300  inches.   The  moment  of  inertia  is  5.2  in 

and  Young's  modulus  of  elasticity  is  30  x  10   psi.   It  should  be 

noted,  however,  that  a  column  of  any  constant  cross  section  or 

material  could  have  been  used. 

Case  I  -  As  illustrated  in  Fig.  lU,  the  pipe  is  subjected  to 
thermal  strain  only.   Curve  1  shows  the  initial  configuration. 
Subsequent  numbered  curves  represent  the  configuration  correspond- 
ing to  increasing  values  of  thermal  strain. 

Case  II  -  Here  the  pipe  is  subjected  only  to  an  increasing 
moment.   Figure  15  shows  the  sequence  of  resulting  configurations. 

Case  III  -  This  example  shows  the  resulting  configuration  of 
the  pipe  (see  Fig.  16)  with  simultaneously  increasing  thermal  strain 
and  left  end  moment. 
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V.   RECOMMENDATIONS 

If  future  studies  are  made  using  the  analysis  set  forth  in  this 
thesis,  the  follovring  areas  are  recommended. 

The  problem  should  be  extended  to  include  multi-span  columns. 
In  the  analysis  of  two  or  more  spans  it  would  be  convenient  to  have 
a  better  way  of  obtaining  the  junction  moment  or  moments.   Perhaps 
a  method  can  be  developed  that  does  not  require  an  initial  approxi- 
mation. 

Another  area  of  study  would  be  to  develop  a  proof  that  snap- 
through  cannot  take  place  or  develop  a  clear  delineation  of  the 
circumstances  under  which  it  could  take  place.   The  study  probably 
should  focus  initially  on  t\-JO-span  columns. 
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APPENDIX  A  -  DESCRIPTION  OF  COMPUTER  PROGRAM 

A.  GENERAL  REMARKS 

The  digital  computer  program  used  to  carry  out  the  analysis  set 
forth  in  this  thesis  is  a  double  precision  FORTRAN  language  program. 
It  is  self  contained  except  for  one  outside  subroutine  used  to  draw 
the  configurations.   This  subroutine  or  a  similar  one  should  be 
available  in  any  computer  library.   Various  write  statements  are 
interspersed  throughout  the  program  to  allow  the  user  to  follov;  the 
progress  of  the  computations,  and  to  help  trace  any  errors.   These 
statements  may  of  course  be  left  out  of  the  program  to  lessen  the 
amount  of  unneeded  output. 

B.  BASIC  FUNCTIONS  PERFORMED  BY  THE  PROGRAM 

The  basic  program  for  the  single-span  plots  the  initial  deflection, 
determines  the  axial  load  on  the  column,  and  plots  the  final  deflection. 
The  axial  load  is  determined  using  a  Newton  Raphson  technique  on  Eq. 
(19).   Iteration  continues  until  convergence  is  reached. 

For  the  two-span  problem,  a  process  is  employed  of  repeatedly 
selecting  values  of  the  junction  moment.   Then  .the  two  spans  are 
analyzed  and  the  two  junction  slopes  are  computed  and  compared.   Then 
the  correct  value  of  the  junction  moment  is  found  using  the  method  of 
successive  decimation.   Although  the  slope  intercept  method  is  more 
accurate  when  it  works,  its  unreliability  causes  it  to  be  less 
desirable  for  use  in  the  general  program.   After  the  junction  moment 
is  found,  the  final  deflection  is  then  plotted. 
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C.   INITIAL  DEFLECTION 

The  Fourier  coefficients  for  the  initial  deflection  of  the  column 

are  computed  for  the  analysis  previously  described.   The  program 

accepts  coefficients  describing  a  polynomial  curve  which  passes 

through  the  left  support  of  span  one  and  the  right  support  of  span 

two.   Then  the  modification  of  this  curve  is  computed  which  enforces 

the  intermediate  constraint;  this  modified  curve  is  also  in  the  form 

of  a  polynomial.   This  yields  two  polynomial  curves,  one  for  each 

span;  cf.  Eqs.  (27)  and  (28).   Then  the  Fourier  coefficients  a  are 

evaluated  for  each  span. 

1.  Subroutine  SINCOE 

This  subroutine  evaluates  the  integral 

IT 

GA(m,n)  =  I  y   sin/nyjdy  (A.l) 

0 
which  is  used  in  the  evaluation  of  the  Fourier  coefficients.   The 

method  involves  evaluating  integral  number  6,  page  117,  of  Ref.  ^  . 

See  program  listing  in  Appendix  C. 

2.  Evaluation  of  Fourier  Coefficients 

The  Fourier  coefficients  are  evaluated  according  to  the 

following  formula: 

L 

a^  =  \    rg(x)  sin  pp^j  dx,  n  =  1,2,...,  (A. 2) 

T) 
where  L  is  the  length  of  the  span  and  n  is  the  number  of  coefficients 

desired.   For  the  left  span,  y,  as  given  by  Eq.  (27),  is  substituted 
for  the  value  of  g(x),  and  for  the  right  span  Eq.  (28)  is  used. 

The  method  carried  out  in  the  program  can  best  be  illustrated 
through  the  use  of  an  example.   Suppose  the  initial  configuration  f(x) 
is  given  by  the  equation 

f(x)  =  Ax^  +  Bx.  (A. 3) 
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Then  Eqs.  (27)  and  (28)  become 
y  =  Ax^  +  Bx  - 


eETOj+L^) 


a^+L2)2  -  x2  -  L22 


] 


CA.H) 


and 


y  =  A(x.UL^)2+  BCx^-^+L^)-  s^;,^,^^^  <=^'-^4)  ['V^z'  "  ^^^'^^l^' 


-  L, 


respectively.   Equation  (A. 4)  can  be  rewritten  as 

y  =  Ax^  +  Bx  -  Cx  +  Dx^ 
and  Eq.  (A. 5)  can  be  rewritten  as 

y  =   A(x''-+L^)^+  B(x-VL^)  -  C(x"+L^)  +  D(x-'--+L  )^  -  Fx-^  , 
where 


and 


C  = 


D  = 


F  = 


WL2(L^+2L^L2) 
GEKL^+L^)    ' 


GEKL^+L^)   ' 


W 


5EI 


For  the  first  span,  Eq.  (A. 2)  becomes 


^n  ^      L 


f     2     .     /niTx\  ,     ^  n   r         •     /nTTxN 

I    x     sm  I  — —  I  dx  +  B   /    X   sm  (  — —  J 

4        ^1^       ■'0       V  1/ 


dx 


The   integrals   are   all  of  the   form 


/m      .     /nTTx  \ 


dx 


(A. 5) 


(A. 6) 


(A. 7) 


(A. 8) 


(A. 9) 


(A. 10) 


(A. 11) 


(A. 12) 
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Making  the  substitution  y  =  -; —  ,  this  integral  may  be  put  in  the  form 

1 

\—WJ         I     y™  sin(nyjdy  =  p""""   GA(m,n)  ,  (A. 13) 

0 
where 

P  =  —  ,  (A. 14) 

and  the  definition,  (A.l),  is  used.   Thus  Eq.  (A. 11)  can  be  rewritten 
in  the  form 

a  =  ^  I  AP^GA(2,n)  +  BP^GA(l,n)  -  CP^GA(l,n)  +  DP^GA(3,n) 


'-'k[ 


=  ^  I   h.    P^"^^  GA(i,n)  (A. 15) 

^1    ^ 

Note,  from  Eqs .  (27)  and  (28),  that  the  maximum  degree  of  the 
adjusted  polynomials  is  Max(3,m"),  where  m^  is  the  degree  of  f(x).   The 
number  n  of  Fourier  coefficients  to  be  employed,  and  the  integer  m'" 
are  among  the  inputs  to  the  program.   Internally,  m""  is  increased,  if 
necessary,  to  a  minimum  value  of  3.   Then  SINCOE  is  used  to  generate 
the  numerical  constants  GA(i,j),  i  =  1,2,  ...,  m''";  j  =   1,2,  ...,  n. 
Thus,  the  coefficients  are  available  for  the  evaluation  of  the  Fourier 
coefficients  as  represented  by  such  equations  as  (A. 15),  the  coef- 
ficients A,  B,  C,  ...,  in  which  are  obtained  in.  an  obvious  manner. 

The  Fourier  coefficients  for  the  second  span  are  found  in 
basically  the  same  manner,  except  that  now  the  computations  are  more 
numerous  due  to  the  necessity  of  expanding  powers  of  (x"+L).   It  can 
be  seen  that  for  the  second  span  Eq.  (A. 2)  becomes  quite  lengthy  since 
g(x)  will  contain  many  terms.   Looking  only  at  the  first  term  in  the 
example,  we  have,  after  expansion 
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'0 
dx 


(A. 16) 


0 
This  can  also  be  written  in  the  form 

a  =  I-  Ja  p2  GA(2,n)  +  L.P^GACl.n)  +  L, ^pGA(0,n)|+  .. A       (A. 17) 


i^   =     I~  ^^  p2  GA(2,n)  +  Lj^P^GAdjn)  +  Lj^2pGA(0,n) 


where  now 

L 


2 
P  =  —  . 


D.  VARIATIONS  IN  THE  PROGRAM 

Different  combinations  of  strain  and  end  moments  to  be  applied  to 
the  column  can  be  incorporated  into  the  program.   This  can  be  accom- 
plished easily  through  the  use  of  the  input  cards  described  below. 
The  program  can  also  be  altered  through  an  additional  DO-loop  to  show 
a  "time  history"  of  the  column.   The  strain  and /or  end  moments  can  be 
incremented  in  the  DO-loop  so  that  intermediate  configurations  may  be 
plotted.   This  is  useful  in  looking  at  snap-through.   The  program 
listed  in  Appendix  C  has  this  variation  incorporated. 

E.  USE  OF  THE, PROGRAM 

This  program  may  be  used  by  inserting  the  follov;ing  21  input  cards 
in  the  designated  places  at  the  beginning  of  the  program.   Each  card 
contains  only  one  entry.   Note  that  these  cards  are  made  a  part  of  the 
program  deck  and  are  not  read  in  as  data  follovring  the  program.   See 
program  listing  for  example.   All  cards  must  be  punched  beginning  in 
column  seven. 


41 


DIMENSION  GA(MPOL,NFOR)  -  integer  numbers  must  be 

inserted  for  MPOL  and  NFOR 

MPOL  =  maximum  degree  of  initial  deflection  function 

NFOR  =  number  of  Fourier  coefficients  desired 

NDEC  =  number  of  decimations  in  determining  junction 
moment 

C6  =  initial  function  coefficient 

C5  =  initial  function  coefficient 

C4  =  initial  function  coefficient 

C3  =  initial  function  coefficient 

C2  =  initial  function  coefficient 

CI  =  initial  function  coefficient 

LI  =  length  of  left  span 

L2  =  length  of  right  span 

E  =  Young's  modulus  of  elasticity 

W  =   cross  section  area  of  the  column  m  inches 

h 
Z  =  moment  of  inertia  of  the  column  m  inches 

EPSl  =  applied  strain  in  left  span  in  inches/inch 

EPS2  =  applied  strain  in  right  span  in  in/in 

SS  =  applied  left  end  moment  in  pound  inches 

UU  =  applied  right  end  moment  in  pound. inches 

XMOMT  =  starting  moment  for  successive  decimation 
iteration 

DELTA  =  initial  increment  for  successive  decimation 
iteration 

Notes : 

1.  All  inputs  are  floating  point  numbers  (D  format)  except  MPOL 
and  NFOR,  which  are  integer  numbers. 

2.  There  must  be  an  input  card  for  each  input  listed,  even  if  it 

is  zero. 
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3.  Each  initial  function  coefficient  must  correspond  to  the 
power  of  X  it  is  associated  with.   That  is 

f(x)  =  C5x^  +  C5x^  +  C^x"*  +  C3x^  +  C2x^  +  Clx 
(A  polynomial  of  degree  greater  than  six  can  be  accommodated  only  by 
making  internal  changes  in  the  program:   these  s?iould  be  rather 
obvious. ) 

4.  The  program  is  dimensioned  for  a  maximum  number  of  50  Fourier 
coefficients  and  a  maximum  combined  span  length  of  400  inches.   (This 
limitation  is  related  to  the  fact  that  deflections  are  computed  at 
two-inch  intervals  and  are  stored  prior  to  output.)   If  the  user  wishes 
to  extend  these  limits,  the  dimension  statements  must  be  altered 
accordingly. 

5.  Likewise  the  call  statement  for  the  DRAV/  subroutine  may  have 
to  be  altered  to  accommodate  the  particular  problem  in  question.   (For 
example,  to  provide  a  desired  scaling.) 

6.  Where  to  start  the  iteration  for  the  method  of  successive 
decimation,  i.e.,  selection  of  the  input  XMOMT ,  poses  a  problem. 
Experience  and  perhaps  a  trial  run  will  best  produce  a  suitable  start- 
ing point. 
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APPENDIX  B  -  FLOW  CHART 


AAAA 


( START  ) 

1 

' 

DIMENSION 

I  DAT A  INPUT 


XCOUNT  =0.0 


ACOUNT  =0.0 
BCOUNT  =0.0 


± 


ACOUNT  =  ACOUNT  +  1.0 


-5»- 


XCOUNT  =  XCOUNT  +  1.0 


1 


COMPUTE  FOURIER 
COEFFICIENTS,  a^ 


-^ 


V 


WWW 
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AAAA 


WWW 


COMPUTE  AND  PLOT 
INITIAL  DEFLECTION 


<- 


V    SI    y 


COMPUTE  LOAD 
COMPUTE  FOURIER  COEFFICIENTS,  b 
COMPUTE  END  SLOPES 


=  1 


KK 


--2 


Mf^ 


=  2 


i 


-A 


COMPUTE  END  SLOPE  ERROR 

USE  METHOD  OF  SUCCESSIVE 

DECIMATION 


BCOUNT  =  BCOUNT+  I .0 
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COMPUTE  AND  PLOT 
FINAL  DEFLECTION 
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